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Abstract 

Making use of the primordially isocurvature fluctuations, which are generated in 
inflationary models with multiple scalar fields, we make a phenomenological model 
that predicts formation of primordial black holes which can account for the massive 
compact halo objects recently observed. 



In inflationary universe models PJ with multiple scalar fields, not only adiabatic 
but also primordially isocurvature fluctuations || are generated in general. For 
example, if we consider ordinary slow roll-over inflation models in Brans-Dicke 
gravity ||, the Brans-Dicke dilaton acquires isocurvature fluctuation out of quan- 
tum fluctuation in the inflationary regime. In this particular model, however, the 
isocurvature fluctuation does not play any important role in structure formation 
because the energy density of the Brans-Dicke field remains as small as ~ oj~ x 
times the total energy density, where u > 500 is the Brans-Dicke parameter || . 

On the other hand, primordially isocurvature fluctuations can be cosmologically 
important if energy density of its carrier becomes significant in a later epoch []3| . 
Furthermore it is relatively easier to imprint a nontrivial feature on the spectral 
shape of the isocurvature fluctuations. Making use of this property here we con- 
struct a model which possesses a peak in the spectrum of total density fluctuation 
at the horizon crossing. Then we can produce a significant amount of primordial 
black holes (PBHs) around the horizon mass scale when the fluctuation at the hori- 
zon crossing becomes maximal. To be specific, we choose the model parameters so 
that one can account for the observed massive compact halo objects (MACHOs) 
in terms of PBHs thus produced. That is, our goal here is to produce PBHs of 
mass ~ O.IMq with the abundance of ~ 20% of the halo mass, corresponding to 
about 10~ 3 times the critical density [§]. 

PBHs are formed if initial density fluctuations grow sufficiently and a high 
density region collapses within its gravitational radius. First let us review its 
formation process. The background spacetime of the early universe dominated by 
radiation is satisfactorily described by the spatially-flat Friedmann universe, 

ds 2 = -dt 2 + a 2 (i) [dr 2 + r 2 (d6 2 + sin 2 Odcp 2 )] , (1) 

whose expansion rate is given by 

with p{t) being the background energy density and a dot denotes time derivation. 
Following Carr ||, let us consider a spherically symmetric high density region 
with its initial radius, R(to), larger than the horizon scale ~ t . The assumption 
of spherical symmetry will be justified below. Then the perturbed region locally 
constitutes a spatially closed Friedmann universe with a metric, 

dr 2 



ds 2 = -dt' 2 + R 2 {t') 



Kr 2 



+ r 2 (d6 2 + sin 2 8d<p 2 ) 



k>0, (3) 
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Then the Einstein equation reads 



, _ , ldR\ 2 8-C; 



H =—"^-mF Y (4) 

there, where p+ is the local energy density. One can choose the coordinate so 
that both the background and the perturbed region have the same expansion rate 
initially at t — t' — to. Then the initial density contrast, 5q, satisfies 



r _ P+0 — PO _ /_x 
PO -Ho-Ko 



where a subscript implies values at t . It can be shown that the two time variables 
are related by [|T0j 

_ , 3 c/t' (it 

The perturbed region will eventually stop expanding at R = _R C , which is 
obtained from 

= ^ P0 (l + *,) f - 4 = i/„ 2 (l 4- *,) ( - (7) 

as 



3 ruv ' UJ \Rj R\ uv ' U/ Vi? c ; Rl 



corresponding to the epoch 



Rc = jii^flo = i„'flo, (8) 

The perturbed region must be larger than the Jeans scale, Rj, in order to contract 
further against the pressure gradient, while it should be smaller than the horizon 
scale to avoid formation of a separate universe. We thus require 

Rj~c s t c < R c < t c , (10) 

or 

Cs < ^ - ^4 < 1, (11) 

where c s is the sound velocity equal to 1/V3 in the radiation dominated era. Since 
RoS /to is time independent, it suffices to calculate the constraint on 5 at a 
specific epoch, say, when the region enters the Hubble radius, 2R = It. We find 
that the amplitude should lie in the range 

\ < 6(R = t) < 1. (12) 
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The gravitational radius, R g} of the perturbed region in the beginning of con- 
traction with R c ~ c s t c is given by 

R g = 2GM~H 2 R 3 c ~^~c 2 s R c < R c . (13) 

This is somewhat smaller than R c but it also implies that a black hole will be 
formed soon after the high density region starts contraction. Thus we expect that 
a black hole with a mass around a horizon mass at t = t c will result and it has 
in fact been shown by numerical calculations (TTTJ that the final mass of the black 



hole is about O(10 ±0 ' 5 ) times the horizon mass at that time. It has been discussed 
that these black holes do not accrete surrounding matter very much and that their 
mass do not increase even one order of magnitude ||. Note also that evaporation 
due to the Hawking radiation is unimportant for M ^> 10 15 g. 
Since the horizon mass at the time t is given by 

M hor = 10 5 (-?-) M , (14) 



lsec, 

what is required in order to produce a significant number of PBHs with mass 
M ~ 0.1M & is the sufficient amplitude of density fluctuations on the horizon scale 
at t ~ 10 _6 sec. Because the initial mass fraction of PBHs, (3, is related with the 
present fraction f^BH as 

P = fl(10 " 6S ' C) ^BH = 10- 8 fiBH, (15) 

where t eq is the equality time, only an extremely tiny fraction of the universe, 
(3 ~ 10 -11 should collapse into black holes. 

That is, the probability of having a density contrast 1/3 < <5< Ion the 
horizon scale at t ~ 10 _6 sec should be equal to (3. Let us assume density fluc- 
tuations on the relevant scale obey the Gaussian statistics with the dispersion 
^bh *C 1, which is the correct property for the model introduced here. Then the 
probability of PBH formation is estimated as 




(16) 
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which implies that we should have 

5bh ^ 0.05, (17) 

to produce appropriate amount of PBHs. Although it is true that an exponential 
accuracy is required on the amplitude of fluctuations in order to produce a desired 
amount of PBHs we have not been able to obtain the correspondence between 5bh 
and f^BH with such an accuracy because numerical coefficients appearing in the 
above expressions, such as 18 in (|T6|), have been calculated based on a rather qual- 
itative discussion. We therefore will not attempt exceedingly quantitative analysis 
in what follows. 

Note also that for 5bh — 0.05 the threshold of PBH formation, 5 = 1/3, cor- 
responds to 6.4 standard deviation. It has been argued by Doroshkevich |12| that 
such a high peak has very likely a spherically symmetric shape. Thus the assump- 
tion of spherical symmetry in the above discussion is justified and it is also expected 
that gravitational wave produced during PBH formation is negligibly small. 

Since the primordial amplitude of density perturbations on large scales probed 
by the anisotropy of the background radiation |13| is known to be 5 ~ 10~ 5 , the 
primordial fluctuations must have such a spectral shape that it has an amplitude 
of 10 -5 on large scales, sharply increases by a factor of 10 4 on the mass scale of 
PBHs, and decreases again on smaller scales at the time of horizon crossing. It is 
difficult to produce such a spectrum of fluctuations in inflationary cosmology with 
a single component. 

In generic inflationary models with a single scalar field 0, which drives inflation 
with a potential V [</>], the root-mean-square amplitude of adiabatic fluctuations 
generated is given by 

on the comoving scale r(<p) when that scale reenters the Hubble radius f2|. The 
right- hand- side is evaluated when the same scale leaves the horizon during inflation. 
Because of the slow variation of <f> and rapid cosmic expansion during inflation, ([18]) 
implies an almost scale-invariant spectrum in general. Nonetheless one could in 
principle obtain various shapes of fluctuation spectra making use of the nontrivial 
dependence of 8 (r(0)) on V[<f>] fl4jl . In order to obtain a desirable spectrum for 
PBH formation with a mountain on a particular scale, we must employ a scalar 
potential with two breaks and a plateau in between [ 15[| . Such a solution is not 
aesthetically appealing. 
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Here we instead consider an inflation model with multiple scalar fields in which 
not only adiabatic but also primordially isocurvature fluctuations are produced. In 
fact it is much easier to imprint nontrivial structure on the isocurvature spectrum 
as mentioned in the beginning. 

We introduce three scalar fields 0i, 02, and 03 in order to generate the desired 
spectrum of density fluctuation. 0! is the inflaton field which induces the new 
inflation [TO with a double-well potential, starting its evolution near the origin 
where its potential is approximated as 

Ufa] =V - ^* (19) 
so that the Hubble parameter during inflation, Hi, is given by 

R 2 _ 8*Vo 

On the other hand, 02 is a long-lived scalar field which induces primordially isocur- 
vature fluctuations that contribute to black hole formation later. Finally 03 is an 
auxiliary field coupled to both 0x and 02 and it changes the effective mass of the 
latter to imprint a specific feature on the spectrum of its initial fluctuations. 
We adopt the following model Lagrangian. 

£ = ^0i) 2 + \(dct>2? + \{dfo) 2 - Vfa, 02, 3 ] + Ant, (20) 

where £ int represents interaction of 0j's with other fields. Here V[0i, 02, 03] is the 
effective scalar potential which is a polynomial of 0/s only up to the fourth order 
except for the inflaton sector £/[0i]: 

Vfa, 02, 3 ] = Ufa] + \{<t>l- 0^)03 + ^03 - ^0302 + ^02 + ^202, (21) 

where Xj, e, u, 0i c , and m 2 are positive constants. m 2 is assumed to be much 
smaller than the scale of inflation, Hi, and it does not affect the dynamics of 02 
during inflation. Hence we ignore it for the moment. 

Let us briefly outline how the system evolves before presenting its detailed 
description. In the early inflationary stage 0i is smaller than lc , and 03 has 
its potential minimum off the origin. Then 02 also settles down to a nontrivial 
minimum, where it can have an effective mass larger than Hi so that its quantum 
fluctuation is suppressed. As 0i becomes larger than lc , 3 rolls down to its 
origin. Then the potential of 02 also becomes convex and its amplitude gradually 
decreases due to its quartic term. However, since its potential is now nearly flat, 
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its motion is extremely slow with its effective mass smaller than Hj well until 
the end of inflation. In this stage quantum fluctuations are generated to 02 with 
a nearly scale-invariant spectrum. Thus the initial spectrum of the isocurvature 
fluctuations has a scale-invariant spectrum with a cut-off on a large scale. 

After inflation, the Hubble parameter starts to decrease in the reheating pro- 
cesses. As it becomes smaller than the effective mass of 2 , the latter starts rapid 
coherent oscillation. 02 dissipates its energy in the same way as radiation in the 
beginning when its oscillation is governed by the quartic term. But later on when 
A 2 02 becomes smaller than m|, its energy density decreases more slowly in the 
same manner as nonrelativistic matter. Thus 02 contributes to the total energy 
density more and more later, which implies that the total density fluctuation due 
to the primordially isocurvature fluctuations or 2 grows with time. Since what is 
relevant for PBH formation is the magnitude of fluctuations at the horizon cross- 
ing, we thus obtain a spectrum with a larger amplitude on larger scale until the 
cut-off scale in the initial spectrum is reached, that is, it has a single peak on the 
mass scale of PBH formation. 

In the above scenario we have assumed that 02 survives until after the PBH 
formation. On the other hand, were 02 stable, it would soon dominate the total 
energy density of the universe in conflict with the successful nucleosynthesis. As 
a natural possibility we assume that 02 decays through gravitational interaction, 
so that it does not leave any unwanted relics with its only trace being the tiny 
amount of PBHs produced. 

Having stated the outline of the scenario, we now proceed to its detailed descrip- 
tion and obtain constraints on the model parameters to produce the right amount 
of PBHs on the right scale. First we consider the evolution of the homogeneous 
part of the fields. During inflation, the evolution of the inflaton is governed by 
the U[4>i] part of the potential. Solving the equation of motion with the slow- roll 
approximation, 



where 0i, is the field amplitude at some initial epoch ^. The above approximate 
solution remains valid until |C/"[0i]| becomes as large as 9H] at t = tf, when 
inflationary expansion is terminated and we find 4>\{tf) = 3H]/Xi. Then (p3|) can 




(22) 



we find 




(23) 
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also be written as 

Al0l(t) = 2H!{t f -t) + \ = 2r(t) + 1 " MTy (24) 

where r(£) is the e-folding number of exponential expansion after t (<tf) and the 
last approximation is valid when r ^> 1. From now on, we often use r(i) as a new 
time variable or to refer to the comoving scale leaving the Hubble radius at t. Note 
that it is a decreasing function of t. 

As stated above, we are taking a view that 0i determines fate of 03 and that 
03 controls evolution of 02 but not vice versa. In order that 03 does not affect 
evolution of 0i the inequality 

AiA 3 > e 2 , (25) 
must be satisfied, while we must have 

A 2 A 3 > v\ (26) 

so that 02 does not affect the motion of 3 . We assume these inequalities hold 
below. 

When 0! < lc , both 2 and 03 have nontrivial minima which we denote by 
2m and 3m , respectively. From 

V 2 = -W> 2 2 + A202 5 = 0, (27) 

and 

V 3 = e(0 2 - 2 c )0 3 + A303 5 - W 2 3 = 0, (28) 
with Vj = dV/d(j)j, we find 

A 2 0LW = W>LW (29) 
A 3 2 m W = e(0 2 c -0 2 (t)) + z/0 2 m (t) = e(0L-0 2 (t)), (30) 

where (|26D was used in the last expression. In the early inflationary stage when 
0i "C 0ic, the effective mass-squared of <f>j, Vjj, at the potential minimum is given 
by 

V2 2 [0i, 2m , 3m ] = 2A 2 2 m = ^(0 2 c - 2 ) ~ ^0 2 c = (31) 

^33 [01 5 02m; 03m] = T Hj, (32) 

AiT c 
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where r c is the epoch when (f) 1 = lc . We choose parameters such that 



ue > AiA 3 t c , and e > \\t c . 



(33) 



Then V22 and V33 are larger than H] initially at the potential minimum, so that 
both 02 and 03 settle down to 02 m (t) and 03 m (£), respectively. 

03m (t) decreases down to zero at r = r c when V3 3 [0 3m ] also vanishes. Then 
^33[03m = 0] starts to increase according to 0i and soon acquires a large positive 
value, which implies that 03 practically traces the evolution of 03 m (t) down to 
zero without delay. On the other hand, 02 evolves somewhat differently because 
it does not acquire a positive effective mass from 3 at the origin. Although 2 (t) 
traces 02 m (t) initially, as V22[02m] becomes smaller it can no longer catch up with 
02m(^)- From a generic property of a scaler field with a small mass in the De Sitter 
background, one can show that this happens when the inequality 



d(j) 2m {t) 



02m (t) dt 



> 



2m 



3H 



gets satisfied, or at 



with 



t= 1 + 



'A1A3 
2z/e 



T c = Ti = T c , 



02 ( T l) = 4>ll = 



ve 3Hf 



(34) 



(35) 



(36) 



2A1A3 2\ 2 n 

Thus 02 slows down its evolution. In the meantime 03 vanishes. Then 02 is 
governed by the quartic potential. We can therefore summarize its evolution during 
inflation as 



02m(^) = aSffc - r) ' T ^ T '' 



<f>\ 



, , 2A 2 ^ 

1 + ~mf 



(ti ~ r) 



(37) 



< T < Tl. 



Since A202; is adequately smaller than H], 02 remains practically constant in the 
latter regime. Let us also write down time dependence of its effective mass-squared 
for later use. 



A1A3 \T C t) 



T > Ti, 



^22 [02] = < 



3A 



2921 



2\ 2< t> 



-1 -1 



(38) 



l + -Sff(7J-T) 



, < T < Tl. 
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Next we consider fluctuations in both the scalar fields and metric variables in a 
consistent manner. We adopt Bardeen's gauge- invariant variables <&a an d &h 103 > 
with which the perturbed metric can be written as 



ds 2 = -(1 + 2<& A )dt 2 + a(t) 2 (l + 2$ H )c/x 2 , (39) 

in the longitudinal gauge. In this gauge scalar field fluctuation 5(j)j coincides with 
the corresponding gauge-invariant variable by itself. 

Assuming an exp(ikx) spatial dependence and working in the Fourier space, 
the perturbed Einstein and scaler field equations are given by 

$a + $h = (40) 
$h + H$ h = -AnGifaSfa + fcSfc + foSfo) (41) 



54>j + SHS'fa + + Vjj 5(f), = 2Vj$ A + $ A fa - 3Qh4>3 - E VjiSfa, (42) 

Note that all the fluctuation variables are functions of k and t. 

The above system is quite complicated at a glance. However, using constraints 
on various model parameters we have obtained so far, i.e. (|25l), (p6|) , and (03|) , 
it can somewhat be simplified. First, since 03 has an effective mass larger than 
H 2 during inflation except in the vicinity of r = r c , quantum fluctuation on 5</> 3 
is suppressed and moreover its energy density practically vanishes by the end of 
inflation. We can therefore neglect fluctuations in 3 . On the other hand, we can 
show that 

■\<Pi\ « N, (43) 




with the help of (p6|) and (0). Hence ([11]) and (fl"2" ) with j = 1 reduce to 



$h + H$h = -4^01(501 (44) 



<$! + 3H5^ + + Vuj 50i = 2V&A - 4$h0i. (45) 

Thus only 0! contributes to adiabatic fluctuations and it can be calculated in the 
same manner as in the new inflation model with a single scalar field. This is as 
expected because 0i dominates the energy density during inflation. In fact since 
we are only interested in the growing mode on the super-horizon regime which 
turns out to be weakly time-dependent as can be seen from the final result, we can 
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consistently neglect time derivatives of metric perturbations and terms with two 
time derivatives in (03) and (H5|) during inflation. We thus find 



H 



4ttG • 



H, 



(46) 



The resultant amplitude of scale-invariant adiabatic fluctuations depends on Ai 
and one can normalize its value using the COBE observation [13[ 



as 



Ax = 1.3 x 10 



-13 



(47) 



On the other hand, 502 satisfies (|42"D with j = 2. From quantum field theory 
in De Sitter spacetime, it has an RMS amplitude 502 — (H 2 /2k 3 ) 1 / 2 when the 
/c-mode leaves the Hubble radius if V22 is not too large. Since V2 vanishes when 
02 = 02m and V2$a remains small even for r < 77, we can neglect all the terms in 
the right-hand side, to yield 



502 + 3^502 + ^ 22 502 = 0, 



(4£ 



when <C a(t)Hi. We can find a WKB solution with the appropriate initial 
condition, 



502 (M) 



2k 3 \S(t)) GW ,n, 



S(t')H 



, - 1*, 



dt'Y 



(49) 



W - 3 1 41/22 

where tk is the time when /c-mode leaves the Hubble radius: k = a(tk)Hj. The 
above expression is valid when \S\ <C S 2 . In terms of r, (f49|) can be expressed as 



50 2 (fc,r) = \ — ( - 7T7 - V ) exp 



S(r) 



2fc 3 V 5(r 



1 



Tk S(r')rfr'-^(r k -r) 



(50) 



I 1 



2A2 0| 



1 + 



<(r* - r) 



1 

-i\ 2 



, < r < n , 



where r k = r(tk) and r* = 3A ^ 3 • The above equality is valid until the end of 
inflation at t = or r = 0. 
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Let us assume the universe is rapidly and efficiently reheated at t = tf for 
simplicity to avoid further complexity, (see, e.g. ]18[ for recent discussion on 
efficient reheating.) Then the reheat temperature is given by 



(51) 



If there is no further significant entropy production later, one can calculate the 
epoch, t(L), when the comoving length scale corresponding to L pc today left the 
Hubble horizon during inflation as 



t(L) = 37 + In 



lpc 



1 



In 



10 10 GeV 



(52) 



Then the comoving horizon scale at t — 10 6 sec, or L = 0.03 pc corresponds to 



r 



= 34 = r m and the present horizon scale ~ 3000Mpc to r = 59. 
On the other hand, fait) and 8 fa (k, t) evolve according to 



fa + 3Hfa + A 2 2 + m 2 fa 
h 2 + 3H5fa + (3A 2 $ + m\)8fa 



0, (53) 
0, (54) 

l/2t, and the latter 



where the Hubble parameter is now time- dependent: H 
equation is valid for k aH. 

When H 2 becomes smaller than A 2 2 m V)i both 2 and 6 fa start rapid 
oscillation around the origin. Using (|50D one can express the amplitude of the 
gauge-invariant comoving fractional density perturbation of fa, A 2 , as 



A 2 = - 

P2 



in the beginning of oscillation. Here 



5 fa 



>2i 



(55) 



r=0 



>2i 



_ 1 ; 2 . ^2 ,4 1 2/2 
P2 = + "J 'fa + 2 m 202 



(56) 



is the energy density of fa. 

Using the virial theorem one can easily show that it decreases in proportion to 

Thus the amplitude of 2 decreases with a _1 (t). 



ml 



2 r^j 



m 



2- 



a" 4 (t) as long as A 2 </> 2 > 

On the other hand, 8 fa has a rapidly oscillating mass term when A 2 <; 
which causes parametric amplification. We have numerically solved equations ( |53"D 
and ( |54J ) with various initial conditions with \ fa\ 3> \8fa\ initially. We have found 



12 



that in all cases the amplitude of <50 2 remains constant as long as m 2 is negligible. 
Thus A 2 increases in proportion to t 1 ^ 2 in this regime and A 2 becomes as large as 



A ^(A^y (5?) 



while the ratio of p 2 to the total energy density, p tot , which is now dominated by 
radiation, remains constant: 

— = **(tT-) 2 - ( 58 ) 
Ptot \M Pl ) y 1 

As A 2 2 becomes smaller than m 2 , 2 and 50 2 come to satisfy the same equation 
of motion, see (|53|) and fl5"4|), and A 2 saturates to the constant value fl57|). At the 
same time p 2 starts to decrease less rapidly than radiation, in proportion to a~ 3 (t). 
Since A 2 contributes to the total comoving density fluctuation by the amplitude 

A = -^A 2 , (59) 

Ptot 

it increases in proportion to a(t) oc t 1 ^ 2 . In the beginning of this stage, we find 
H 2 = m 4 2 /X 2 (j) 2 { , to yield 



_• 



A-8f^ ^ K=T- . (60) 



02f V m 2 / V^Pi/ \VA20 



2f , 



at a later time £. 

In order to relate it with the initial condition required for PBH formation, we 
must estimate it at the time fc-mode reenters the Hubble radius, t* k) defined by 

k = 2na(tl)H(tl) = (61) 



f) 2 



Since k can also be expressed as k = 27ra/e Tk Hi, the amplitude of comoving 
density fluctuation &t t = t* k is given by 



1 ' , \ 2 



>2i I V A 2 2 f \ 2 / ' '21 



^^^J [wj e " (62) 

In the beginning of this article we discussed the necessary condition on the 
amplitude of fluctuations for PBH formation using the uniform Hubble constant 
gauge. Hence we should calculate the predicted amplitude in this gauge, which is a 
linear combination of A and the gauge-invariant velocity perturbation. However, in 
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the present case in which A grows in proportion to a(t) in the radiation-dominant 
universe, one finds that the latter quantity vanishes and that density fluctuation in 
the uniform Hubble constant gauge coincides with A. Thus we finally obtain the 
quantity to be compared with <5bh in (|ISD, namely, the root-mean-square amplitude 
of density fluctuation on scale r = 2ir/k at the horizon crossing, S(r), as 



5(r) 



Ank 3 
(2n) 3 



\A(k,t 



e Tk C f (r k , r c ,n), 



with 



Cfa, r c , r.) = (j^j exp {£ S(r')dr' - 3 -r k 



where we have used (|50|). We also find 



'3t* 



/ 2 
3t* 



-i 



'3r* 



3Hf 
2A 2 r c ' 



(63) 



(64) 



(65) 



from (|3|) and ([37]). 

The remaining task is to choose values of parameters so that S(r) has a peak on 
the comoving horizon scale at t = 10 _6 sec, which we denote by r m , corresponding 
to Tk = r m = 34, with its amplitude S(r m ) = 0.05. We thus require 



d\nS(r) 



+ S(Tk) " 2 



(66) 



4r k [r c r k + r* (r c - r k )] ' 2 

vanishes at r k = r m , which gives us a relation between r c and r*. Since r c roughly 
corresponds to the comoving scale where scale-invariance of primordial fluctuation 
A 2 is broken, the peak at 5(r m ) becomes the sharper, the closer r c approaches r m . 
For example, if we take r c = 30 we find 

Ave 



so that 



3Ai A 2 



C f = 0.13 and 



200, 



0.045- 



In order to have S(r m ) = 0.05, we require 



% VM, 



'pi 



1.7 x 10 



-15 



(67) 



(68) 



(69) 
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which can easily be satisfied with some reasonable choices of A2 and Hj. However, 
it is not the final constraint. Since we are assuming that the universe is dominated 
by radiation at this time, we require 

92 2. (^l) 2 ( ™} \ ^ e- « 1. (70) 



(hat \MpiJ \V^202f-ff/ / 

Furthermore $2 should decay some time after t = 10 _6 sec so as not to dominate the 
energy density of the universe which would hamper the primordial nucleosynthesis. 
Assuming that it decays only through gravitational interaction, its life time is given 
by 

M 2 Pl 1n _ 5 . B / m 2 \~ 3 



Now we have displayed all the necessary equalities and inequalities the model 
parameters should satisfy. Since there is a wide range of allowed region in the multi- 
dimensional space of parameters, we do not work out the details of the constraints 
but simply give one example of their values with which all the requirements are 
satisfied. 



Hi 


= 1.7 x 10 10 GeV, 


m 2 


= 3.2 x 10 6 GeV, 


Ax 


= 1.3 x 10~ 13 , 


A 2 


= 1.4 x 10" 6 , 


A 3 


= v = 6.7 x 10" 8 , 


e 


= 2.0 x 10~ u , 


= 10" 


6 sec and inequalities 



(72) 



for which p 2 / 'ptot = 0.1 at t — 10 sec and inequalities ( pS] ) and ( p6| ) are maximally 
satisfied. 

Thus we have reached a model with the desired feature making use of a simple 
polynomial potential (pl|). In order to set the order of magnitude of the mass 
scale of the black holes and that of their abundance correctly, we must tune some 
combinations of model parameters such as (|6?f ) and (|69| ) with two digits' accuracy. 
However, there exists a wide range of allowed region in the parameter space to 
realize it. 
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